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Abstract 



A universal phase diagram for type-II superconductors with weak point 
pinning disorder is proposed. In this phase diagram, two thermodynamic 
phase transitions generically separate a "Bragg glass" from the disordered liq- 
uid. Translational correlations in the intervening "multi-domain glass" phase 
are argued to exhibit a significant degree of short-range order. This phase 
diagram differs significantly from the currently accepted one but provides a 
more accurate description of experimental data on high and low-Tc materials, 

simulations and current theoretical understanding. 
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Quenched randomness destabilizes the Abrikosov flux-hne lattice in a pure type-II super- 
conductor, yielding novel glassy phases In the Bragg glass (BrG) phase 0, translational 
correlations decay asymptotically as power laws. On increasing disorder or temperature, 
the BrG phase is unstable to a phase in which such correlations decay exponentially [Q. At 
large values of the applied field if, the disordered liquid (DL) evolves smoothly into a glassy 
state Pi when the temperature T is reduced. This glassy state may be a new thermody- 
namic phase, the "vortex glass", with spin-glass-type correlations in the superconducting 
order parameter, separated from the disordered liquid by a line of phase transitions 0]. An 
alternative picture, attractive in many respects, is that this state resembles a structural 
glass 0, topologically disordered at the largest length scales but lacking long-ranged phase 
correlations 0. 

Fig. 1 summarizes the currently popular view of phase behaviour in the disordered mixed 
phase of type-II superconductors The BrG phase is shown to be unstable to the 

vortex glass on increasing if; it also melts directly into the DL phase upon increasing T. 
The continuous vortex-glass to DL transition line meets the BrG-DL transition line at a 
multicritical point 0. Recent studies indicate that the underlying field-driven Bragg 



glass- vortex glass transition is first-order. 

While this phase diagram hy construction refiects early experiments on YBCO and 
BSCCO [H, many systems, among them the low-Tc superconductors 2H-NbSe2 |TT|-[T3|, 
CeRu2 |ll3|-[15| and Ca3Rh4Sni3 [|^, the cuprate NCCO ll^, the bismuthate (K,Ba)Bi03 



T8| , the borocarbide YNi2B2C |T5| and the mercury and thallium pO[ based compounds. 



behave qualitatively differently. Even for the high-Tc systems, more recent experiments pT[ 



suggest features inconsistent with Fig. 1. That these differences imply qualitively different 
types of phase behaviour is one possibility. Alternatively, such differences might simply 
refiect different limits of an underlying universal phase diagram; this possibility is both 
theoretically attractive and economical. It would be sensible, however, to require that Fig. 
1 be recovered from such a phase diagram in the appropriate limit. 

This Letter proposes Fig. 2 as a universal phase diagram for all weakly disordered type-II 



superconductors. The term "multi-domain glass" (MG) describes the narrow sliver of glassy 
phase which we propose always intervenes between BrG and DL phases. This sliver expands 
into the putative "vortex glass" phase of Fig. 1, as if is increased. Experiments indicate 
that the MG phase melts via a first-order phase transition on T scans at intermediate 
H; this first-order line is shown to meet a line of continuous transitions at a tricritical 
point, although more complex topologies are possible. The inset shows the expected phase 
behaviour at if ~ Hd- The predicted ||2^ reentrance of the glassy phase at low fields 



was established in Ref. Recent experiments [T^Jl^ clearly show the disordered phase 



enveloping a quasi-lattice phase at small H. Fig. 2 incorporates recent proposals for such 



phase behaviour in low-T^ superconductors [13|, derived from an analysis of peak effect 



phenomena in these systems [|T5|, the anomalous and sharp increase in critical currents jc, 
seen over a narrow regime below Hc2, the "peak regime", as H ot T is varied. 

The "multi-domain" glass is argued to be an equilibrium phase, distinct from the vortex 
glass, consistent with recent theoretical 0, experimental and simulational understanding. 
Theoretically, screening destroys vortex-glass-type order in three dimensions I^Sj; recent 



experiments |^ question earlier reports of vortex-glass scaling. Simulations of classical 
directed lines interacting with quenched point impurities see topologically disordered glassy 
phases in equilibrium in a regime which intervenes between relatively ordered (BrG) and 



disordered (DL) phases ||2^. Arguments for the equilibrium nature of this phase and a brief 



survey of its glassy attributes in the context of low-T^ materials were provided in Ref. [jT3| 



What symmetry might be broken across the MG-DL transition? A recent study of 
disordered hard-sphere fluids obtains, within mean-field theory, a one-step replica symmetry 
breaking glass transition for sufficiently strong disorder; these authors stress the resemblance 



of the features of the mean-field phase diagram they obtain to Fig. 1 ||28|. Should these 
ideas be applicable to the disordered mixed phase, a discontinuous freezing of the disordered 
density configuration characterizing the glassy phase would be expected, as against relatively 
smooth behaviour of the entropy and internal energy across this transition. An alternative 
possibihty for the MG phase, consistent with some of the proposals here and testable via 



an analysis of Bitter decoration patterns, is a hexatic glass ||2^ ; the recent discovery of 



an equilibrium hexatic phase in a dense system of long oriented DNA molecules favours a 
similar possibility for line vortices. 

We argue that translational correlations in the MG phase can be fairly long-ranged (un- 
like correlations in typical glassy phases), and propose that structure in the MG phase just 
above the BrG-MG phase boundary is best described as a mosaic of ordered domains, with 
typical scale ^ a, the inter-line spacing, for weak pinning. should be largest in 
the interaction-dominated regime but should decrease rapidly for much larger or smaller 
H, reflecting the increased importance of disorder both at high-field and low-field ends 
T3| , p!5|p3| , |3T|| . Some diagnostics for the multi-domain character of this phase are the follow- 



ing: The transition between BrG and DL phases can occur via several intermediate stages, 
as individual domains melt as a consequence of disorder-induced inhomogeneities in the 



melting temperature |^ . Fluctuations of a domain-like arrangement can yield noise signals 
of large amplitude associated with relatively few fiuctuators intermediate structure 

in the ac susceptibility as a function of T jr5| , P^ , a stepwise expulsion of vortices on cool- 
ing and a host of related features [T^,|5S| of the experimental data which are hard to 
rationalize in other pictures. 

The MG-DL transition is interpreted as the true remnant of the underlying freezing tran- 
sition in the pure system. The presence of a critical point on this transition line [§] follows 
from the observation that signals of melting should cease to be seen once Rd is comparable 
to correlation lengths in the pure liquid at freezing. This phenomenology also rationalizes 



O] a large body of data on anomalies associated with the peak effect phenomenon [11 



including "fracturing" of the fiux-line array [|T2[, the association of thermodynamic melting 



with the transition into the disordered liquid |^ as well as conjectures regarding the dy- 
namic coexistence of ordered and disordered phases in transport measurements in the peak 
regime |^ . 



A simple conjecture for the scale of such domain sizes suggests R^ ~ Ra, a Larkin 
length, at least at intermediate H. This length is computed using the coarse-grained free 



energy Fei = f dr±dz [^(V± ■ u)^ + x u)^ + £f(d,u)% where u{r±,z) is the 

displacement field at location (r_L, z) and the integral represents the elastic cost of distortions 
from the ideal crystalline state with vortex lines centred at Rj, governed by the values of the 
elastic constants for shear (cge), bulk (cn) and tilt (C44) Quenched random pinning is 
incorporated by adding Fa = J dv±dz Vd{'r±, z)p{'r±, z), to Fg/, where p{r±, z) = Y,i 5^'^\'r± — 
R— u(Rj, 2;)). This term models the interaction with a quenched Gaussian disorder potential 
Vrf, assumed to be correlated over ^, the coherence length i.e. [Vd(x)Vd(x')] = K{x — x'), 
with K[x — x') a function of range ^ ^j. The notation x denotes (r_L, z). 

The correlator B{r^,z) = [{u{r±,z) — u(0,0))^] where (■) and [■] denote thermal and 
disorder averages respectively, has the following properties 0: For r±,z <C Rc,Ll (the 
transverse and longitudinal Larkin pinning lengths), correlations behave as in the Larkin 
"random force" model i.e. B{x) ~ x'^~'^. At length scales between R^. and the Larkin 
length scale referred to above at which disorder and thermal fiuctuation-induced positional 
fiuctuations become of order a, B{x) ~ x^^™ with Crm ~ (4 — (i)/6 ~ 1/6. The Larkin length 
scale can then be estimated via Ra — RdjY^''™. At still larger length scales, B{x) ~ log(x) 



We estimate the Larkin pinning lengths Rc and from = Rc = 



2 3/2 1/2 



Up 

and Vc ~ R^Lc. Equating BjVc to the energy gain from random pinning f{nVcY^^ yields 
the standard weak-pinning expression jc = -^f (^y-) ^ ■ -^c, > A, jc ~ jdp and 

~ ^Rc H, with jdp the depairing current. For 2H-NbSe2: T^ ~ 7i^, A ~ 700A {H \\ c), 
^ ~ 70 A and a ~ 450A at i? ~ IT. Using id jdp ~ 10~^, we obtain transverse and 
longitudinal Larkin pinning lengths ~ 15 and ~ 24. 

At fields ~ IT, (|)i/C™ _ (^)6 ^ 7 X 10^ leading to Rd/a ~ Ra/a < 10^. Assuming, 
conservatively, Rd ~ lO^a and a longitudinal Larkin length comparable to the size of the 
sample in the c-axis direction, an estimate for the number of "independent fiuctuators" 



obtained in noise measurements on 2H-NbSe2 |3^|T3 can be obtained. For a sample in the 
shape of a thin platelet of transverse area A ~ 1mm x 1mm and with Ad ~ R\, the number 
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of such fluctuators Nj ~ ^ ~ 10^ — 10^, numbers small enough to yield strong non-Gaussian 
effects [p3l . 

Experimentally, behaviour in the mixed phase of low-T^ materials such as 2H-NbSe2, 
CeRu2, Ca3Rh4Sni3, YNi2B2C and several related compounds exhibit a remarkable com- 
monality [|1^,0. In these relatively pure materials [ic/Up ~ 10"'^) where thermal fluctua- 
tions are substantial (Ginzburg numbers Gj ~ 10~^ here, as against typical values for low-Tc 
systems of about 10"®), there is now considerable evidence for a two-step transition from a 
relatively ordered low-T phase into a highly disordered fluid phase This transition 



occurs via an intermediate and profoundly anomalous regime with glassy properties |113 



These experiments indicate that the two transition lines which separate the BrG from the 



DL phase can always be separately resolved. The data on NCCO and (K,Ba)Bi03 118 
are consistent with Fig. 2 but not with Fig. 1. 

Hall-probe based magnetization measurements on BSCCO favour a single-step tran- 
sition into the disordered liquid out of the BrG phase [^]. However, recent muon-spin 



rotation experiments see two transitions ||2T[. The asymmetry of the field distribution, 
a = {AB^)^^^/ (Ai?^)^/^, exhibits a first transition in which a jumps from a value of about 
1.2 (consistent with a vortex-line crystal) to a value of about 1 (indicating a considerable 
degree of local translational order) as T is increased. Across a second transition boundary, 
obtained on further increasing T, a drops abruptly to values close to zero, characteristic of 
the fluid. 

Early experiments on YBCO found that melting in relatively clean untwinned crystals 
was best described as a "complex, two-stage phenomenon" Recent simultaneous mea- 
surements of ac susceptibility and magnetization in this material see a sharp peak effect 
in which the location of the peak in ac susceptibility correlates exactly to the magnetiza- 
tion jump which signals melting. The regime in which x' shows non-trivial signatures of 
the transition is enlarged, compared to the magnetization jump which occurs at a sharply 
defined temperature, a feature which follows transparently from the discussion here - the 
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width of the transition region is related to the width of the shver phase, while the most 
prominent signatures of melting should generically be obtained across a line, the MG-DL 
transition line. 

When does the phase diagram of Fig. 1 approximate that of Fig 2? The sliver of MG 
phase is exceedingly narrow at intermediate H and for weak disorder, a consequence of the 



reduction of effective disorder in an interaction-dominated regime ||T^. Thermal smoothen- 
ing of disorder should reduce the width of this sliver further; thermal fluctuations smear the 
disorder potential Up{T) seen by a vortex line, strongly renormalizing the effective disorder if 
[< >] > S,'^. For a single line pinned by quenched disorder Up{T) ~ f/p(0) exp[— c(T/Trfp)^] 



||7|j2^ where Up(0) is the pinning potential per unit length at T = 0, c is a numerical constant 
and Trfp is the depinning temperature; T^p ~ (Up^^^o/l'^y^^ where cq = ($o/4vrA)^, and 7 
is the mass anisotropy. For YBCO, T^p ~ 20 — 30K |^^. Such a substantial (exponential 



in the simple limit above) reduction in effective disorder should render the sliver unobserv- 
able in many types of experiments, yielding an apparent single melting transition out of the 
ordered phase. 

The ratio R = T^^'^T^'sh/Thighrpiow 

measures the relative importance of thermal fluctu- 
ations in low and high materials; T^ is the melting temperature in the pure case and 
we have approximated T^ ~Tc. Using the following values: NbSe2 '■ Up = 10K/A,C, = 
70A, X = 700A, 7 = 5 and YBCO : Up = lOK/A, ^ = 20A, A = 1400A, 7 = 50 we obtain 
R ~ 10^, suggestive of the relative importance of this effect to the high-T^ materials vis. a 
vis. its irrelevance in low-Tc systems except very close to Hc2. Experimentally, the nomi- 
nal irreversibility line for ~ 3 — 7T in YBCO lies well below the melting line in weakly 
disordered samples Thus, thermal melting occurs in a nearly reversible regime. Local 



Hall probe-based susceptibility measurements on YBCO |^T| find that jc is actually finite 
below Tm but extremely small (~ O.AA/cm'^). These experiments see a peak effect close to 
and below Tm{H), a feature inaccessible in usual SQUID based magnetization experiments. 
Other experiments at somewhat smaller H see a very clear and sharp peak effect with 



properties similar to those in the low-Tc materials ||13|. These results illustrate that sig 



nals of a two-step transition in high-Tc materials may be very hard to access, particularly 
if discontinuities in jc or magnetization across the first transition are small. In contrast, 



experiments on the low-Tc materials discussed in Refs. [|T5|,|T3| indicate an irreversibility line 
located in the fluid phase and provide clear evidence of a two-step transition. 

Anomalously small values of jc translate to anomalously large values of Rd- The 
magnetization discontinuity across the BrG-MG phase boundary should scale roughly as 
the density of unbound dislocations pd in the MG phase for small pd ~ ^/Rl, yielding 
AM ~ AMo{a/Raf where AMq is the magnetization jump in the pure system at T^. Even 
if Ra/a ~ 30, the corresponding induction jump AB ~ AM is of order 10"'^ Mq or within 



noise levels in a typical experiment |]37| , |32| . Thus, for sufficiently weak disorder, the sample 
may behave substantially as a single domain in the intervening MG phase, implying the 
absence of signals of the BrG-MG transition in a magnetization experiment and thus the 
phenomenology of Fig. 1. 

The apparent vanishing of the sliver rationalizes the putative multicritical point 
of Fig. 1. While the sliver may not be resolvable in some classes of experiments, it may 
be apparent in others, particularly those which probe local correlations. However, in the 
thermodynamic limit, provided other, unrelated direct instabilities to the liquid do not 
intervene, continuity suggests two transitions as in Fig. 2. In essence, the phenomenology 
of Fig. 1 assumes by fiat such an instability. There appears to be little justification for 
this assumption; we suggest that it need not hold. For more disordered superconductors, 
whether high-T^ or low-T^, a two-step transition should generically be seen. 



Many recent simulations, for example Ref. |27], see an intermediate field MG phase 



with translational correlations comparable to system sizes, although the dynamics changes 



abruptly across the BrG-MG transition. Interestingly, Ref. ||2^ concludes that the existence 
of a sliver of MG phase always preempting a direct BrG-DL transition cannot be ruled out, 
as in the recent simulations of Sugano et. al. p8[ . 

In conclusion, it is suggested that the BrG phase in all disordered type-II superconductors 
generically transforms first into an intermediate glassy state on heating, rather than directly 
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into a liquid; the proposal of a generic two-step transition out of the disordered liquid 
state should be experimentally testable in thermodynamic measurements. More details will 
appear elsewhere p5| . 

I thank C. Dasgupta for a critical reading of the manuscript. These ideas evolved out of 
many discussions with S. Bhattacharya; his insights were crucial to this work. D.R. Nelson 
is thanked for a useful conversation regarding hexatic glasses. Support from a DST (India) 
Fast Track Fellowship for Young Scientists is gratefully acknowledged. 
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FIGURES 

FIG. 1. The current view of the phase diagram of disordered type-II superconductors 0,^-^. 
In addition to Meissner and normal phases, this phase diagram subdivides the mixed phase into 
three phases - the Bragg glass (BrG), the vortex glass (VG) and the disordered liquid (DL); 
these are described in the text. Note that the Bragg Glass phase melts directly into the liquid at 
intermediate field values. 

FIG. 2. Proposed universal phase diagram for disordered type-II superconductors. The Meiss- 
ner, normal, BrG and DL phases are as shown in Fig. 1. The MG phase (see text) intrudes 
between BrG and DL phases everywhere in the phase diagram. The inset expands the boxed re- 
gion of the main figure, illustrating the reentrance of the BrG-MG phase boundary. Solid lines 
indicate discontinuous transitions. 
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